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In this Letter a A'^-body theory for the radiative heat exchange in thermally non equilibrated dis- 
crete systems of finite size objects is presented. We report strong exaltation effects of heat flux which 
can be explained only by taking into account the presence of many body interactions. Our theory 
extends the standard Polder and van Hove stochastic formalism used to evaluate heat exchanges 
between two objects isolated from their environment to a collection of objects in mutual interac- 
tion. It gives a natural theoretical framework to investigate the photon heat transport properties of 
complex systems at mesoscopic scale. 
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The photon heat tunneling between two bodies has at- 
tracted much attention in the last decades since it has 
been predicted that the heat flux (HF) can exceed, at 
nanoscale, the far field limit set by Planck's black body 
law by several orders of magnitude [l|, y|. This discov- 
ery has opened the way to promising technologies for 
energy conversion and data storage as for example the 
near- field thermophotovoltaics [3|, \^ and the plasmon as- 
sisted nanophotolitography [5|. This dramatic increase 
is generally speaking due to the contribution of evanes- 
cent modes, which are not accounted for in the Stefan- 
Boltzmann law and become only important if the dis- 
tance between the objects is smaller than the thermal 
wavelength [Gi]. The detailed mechanisms which lead to 
such an enhancement are nowadays for a number of ge- 
ometries and materials well understood |6l422| and recent 
experiments [23-27] have confirmed all theoretical predic- 
tions both qualitatively and quantitatively. 

However, some questions of fundamental importance 
remain unsolved in complex mesoscopic systems. Indeed, 
so far, only the HF between two objects [6l-t9| out of 
equillibrium has been considered, but how does the heat 
transport for a collection of individual objects in mutual 
interaction look like? The collective effects in such many 
particle systems has not been explored yet, although it is 
of prime importance for understanding the different heat 
propagation regimes in disordered systems, determining 
the thermal percolation tresholds in random nanocom- 
posites structures and studying thermal effects due to 
the presence of localized modes in such systems. 

Inside a discrete system of bodies maintained at differ- 
ent temperatures the local thermal fiuctuations give rise 
to oscillations of partial charges which, in turn, radiate 
their own time dependent electric field in the surrounding 
medium. These thermally generated fields interact with 
the nearby bodies and modify through different cross in- 
teractions all these primary fields to generate secondary 



fields which in turn affect the radiated fields and so on. 
Generally speaking, this problem belongs to the vast cat- 
egory of many-body problems which constitute the theo- 
retical framework of numerous branches of physics (celes- 
tial mechanics, condensed matter physics, atomic physics, 
quantum chemistry) . A general theoretical framework to 
treat the many-body problem of non-radiative photon 
heat transport does not yet exist. In this Letter we in- 
troduce a self-consistent theory to describe heat transfers 
inside thermally non equilibrated discrete systems. After 
deriving the HF exchanged between two individual ob- 
jects in mutual interaction inside a N-body system, we 
investigate the thermal conductance between a couple of 
particles versus the position of a third object inside a 
three body system. We highlight some emergent phe- 
nomena which specifically result from many-body inter- 
actions. In addition, we will show that for systems with 
at least three objects at different temperatures one can 
actively control the heat flow in nanoscale junctions, i.e., 
one has a thermal heat transfer transistor. 

To start, let us consider a discrete set of N objects 
located at positions r^ and maintained at different tem- 
peratures Ti with i — I, . . . , N. Suppose that the size of 
these objects is small enough compared with the small- 
est thermal wavelength Xt^ — ch/{kBTi) so that all 
individual objects can be modeled to simple radiating 
electrical dipoles. For metals one has also to include 
the magnetic dipole moments due to the induction of 
eddy currents |lO| . Such an extension of our approach 
is straightforward so that for convenience we will con- 
sider electric dipoles only. The Fourier component of 
the electric field at the frequency w (with the conven- 
tion f{t) = / ^ f {^)^~^^*) generated at the position r.i 
by the fiuctuating part p^"'^ of electric dipole moment of 
the particle j which is located at rj reads 
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with Ho the vacuum permeabihty and (G'J = (G(ri,rj;a;) 
the dyadic Green tensor (i.e. the propagator) between 
the particles i and j inside the set of N particles. On the 
other hand, by suming the contribution of fields radiated 
by each particle, the dipolar moment induced by the total 
field on the z-th particle is given by 



Now, for calculating the Green's function (GF) for a 
system of A^ particles we use the set of 3A^ self-consistent 
equations 
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where a^ is the particle's polarizability and Eq is the vac- 
uum permittivity. Then, the power dissipated inside the 
particle i at a given frequency uj by the fluctuating field 
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vacuum GF defined with the unit vector r 
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Vij being the vector linking the center of dipoles i and j, 
while rij —\ Yij |and 1 stands for the unit dyadic tensor. 



Eij generated by the particle j can be calculated from Inserting relation ([T]) into this sytem leads to the GF 



the work of the fiuctuating electromagnetic field on the 
charge carriers as 



Vj-,,,u: = 2Re(-zL^pr(w) • E* 



(3) 



where the brackets represent the ensemble average. Us- 
ing relations ([1]) and ^ between the dipole moments 
and the fiuctuation dissipation theorem, i.e. {p*j aPi,ii) = 
2— Im(aj)8(a;,rj)(5a^5ij we find after a straightforward 
calculation that 

/■OO 1 

introducing the transmission coefficient (TC) 
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T^,3{^) = -— Im(a,)Im(a,)Tr[G''^G''^t]. 
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In order to present the HF in an obvious Landauer-like 
manner [29|,|30| we rewrite the HF in terms of the conduc- 
tance Gij = dPj^,/dTj so that Vj^., ^ G.^AT. Then 
we find 



r,^. = 3 
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, ,, = Jdx f[x)%,j (cc)/(7rV3) is the mean TC [2^ 



where Ti. 

with f{x) = cc^ exp(— x)/(exp(x) — 1)^. This expres- 
sion generalizes the Meir-Wingreen Landauer type for- 
mula l3l| for photon HF in N-body systems. In the 
case of two particles {N = 2) one can easily show that 
7i J (w, d) G [0, 1] and therefore Ti^j £ [0,1] as well. Hence 
the conductance between two dipoles is limited by 3 times 
the quantum of thermal conductance n^kBT/{3h) [32|- 
In other words, only three channels contribute to the HF 
between two dipoles namely the channels due to the cou- 
pling of the three components pj^a with the same three 
components Pj^Q (i.e. same polarization). Of course, by 
adding further particles this limit cannot be exceeded, 
whereas the HF can be increased or decreased with re- 
spect to the case of two particles. Nevertheless, the num- 
ber of channels increases if electric multipoles [33, 34] as 
well as the magnetic moments J10l| come into play. 
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(9) 
With these relations and Eq. ^ at hand it is possible to 
determine the interparticle HF in a system of TV particles 
out of equillibrium. 

Let us now apply this theoretical formalism to describe 
some emerging many-body effects. To this end, we con- 
sider the simplest possible configuration where such ef- 
fects occur that is a triplet of particles. We consider only 
the interparticle HF between particle 1 and 2 separated 
by a distance 21 in the presence of the third particle. 
Here, we assume that Ti = 300 K and T2 = T3 = 0. The 
interparticle HF is then given by 



f 12(21, rs) = 7'i^2 -^2^1 = 7'i^2- 
In this case, the dyadic GF reads 
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^Gl'aiGl^a2, and B^i^ = G^Sag + ^QfaiQ}^'a3. 
Some numerical results are shown in Fig. 1 and Fig. 2. 
We plot the resulting interparticle HF (/)*2 between par- 
ticle 1 and 2 in the presence of body 3 normalized to the 
HF for two isolated dipoles. In both Figs, the position of 
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Figure 1: Normalized HF exchanged between two SiC spher- 
ical particles maintained at Ti = 300K (left particle) and 
r2 = OK (right particle) with respect to the position rs — 
[xz,y-i) of a third SiC particle of same radius and Ta = OK 
for (a) 21 = 600 nm, and (b) 2/ = 800 nm. The HF is normal- 
ized by the HF exchanged between two isolated dipoles in the 
same thermal conditions. The dark zone with a negative HF 
corresponds to the region which cannot be occupied by the 
third particle. For the sake of clarity we consider here that 
all particles are identical (100 nm radius) and their electric 
polarizability given by the simple Clausius-Mossotti form [28] 
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R denoting the particles radius . The dielec- 



tric permittivity of particle is described by a Drude-Lorentz 
model. 



the particles for which the interparticle HF is calculated 
is fixed, but the position of the third particle is changed. 
It can be seen that for some geometric configurations the 
HF mediated by the presence of the third particle can be 
larger than the value we usually measure for two isolated 
dipoles. In particular, we observe an exaltation of HF 
of about one order of magnitude when the third particle 
is located between the two other particles, i.e., when all 
three particles are aligned. Hence, the HF between two 
dipoles can dramatically be increased when inserting a 
third particle in between. Note, that the three body sys- 
tem described above represent a photon heat transistor 
where the heat HF between two particles can be actively 
controlled by the presence of a third particle. This could 
be achieved through classical AFM manipulation tech- 
niques. 

In order to understand some of the physics behind 
this enhancement mechanism, we show in Fig. 3 the TC 
72,1 (t^) for two isolated particles and three aligned par- 




Figure 2: Normalized HF between two spherical particles of 
same radius with respect to the position of a third one which 
is equidistant to both particles (Ti = iQQK , T2 = T3 = OK 



tides for a frequency range around the surface phonon 
resonance of the particles and for different interparti- 
cle distances. First of all, it can be observed that the 
TC shows different resonances where it is close to 1. 
These resonances can be found from the expression for 
the GF of the N-body system in Eq. (|S]) by evaluating 
det(l— ^0) — where Ao is defined in Eq. ^. In fact, for 
the considered systems of two or three aligned particles 
one yields three configurational resonances [35|- Two of 
these resonances are degenerate because of the rotational 
invariance around the alignement axis (for more details 
see Ref. [37|). Apart from these configurational reso- 
nances we have the surface mode resonance at w = ujsr 
with e{u;sr) = —2, which becomes dominant for large dis- 
tances so that one sees only one resonance in this case. 
On the other hand, for distances close to the particle 
radius the multiple interactions become dominant and 
several resonances show up in the TC. Note, that the 
dipole model [71 18| is only valid for I > 2R. Now, from 
Fig. 3 it is clear that at long separation distances the 
coupling between two dipoles becomes more efficient in 
presence of a third mediator than without, so that the 
HF enhancement can be attributed to a three body ef- 
fect that is a resonant surface mode coupling mediated 
by the third particle. Nevertheless, the absolute value for 
the interparticle HF is still far away from the theoretical 
upper limit. 

In conclusion, we have introduced a theoretical frame- 
work to investigate photon heat transport in mesoscopic 
systems where strong electromagnetic interactions exist. 
In particular, a Meir-Wingreen-Landauer-type formula 
for the radiative HF through N-body interacting photon 
regions has been derived. A detailed study of three body 
systems has allowed to identify a many body exaltation 
mechanism of HF due to configurational resonances. This 
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Figure 3: TC between two spherical SiC particles separated 
by a distance 21 (top) and between the same particles when 
a third SiC particle is located on the mass center of this cou- 
ple (bottom). All particles have the same radius. The vertical 
dashed lines marks the distance where the particles are touch- 
ingHi. 



effect could be used to improve for example the perfor- 
mance of near- field thermophotovoltaic devices [3|, |4| , by 
placing nanoparticles on the surface of photovoltaic cells. 
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